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Abstract. We report the fusion of photons from two independent photonic crystal 
fiber sources into polarization entangled states using a fiber-based polarizing beam 
splitter. We achieve fidelities of up to = 0.74 ± 0.01 with respect to the 
maximally entangled Bell state using a low pump power of 5.3mW with a success 
rate of 3.2 four-fold detections per second. By increasing the pump power we find that 
success rates of up to 111.6 four-folds per second can be achieved, with entanglement 
still present in the fused state. We characterize the fusion operation by providing 
a full quantum process reconstruction. Here a model is developed to describe the 
generation of entanglement, including the main causes of imperfection, and we show 
that this model fits well with the experimental results. Our work shows how non-ideal 
settings limit the success of the fusion, providing useful information about the practical 
requirements for an operation that may be used to build large entangled states in bulk 
and on-chip quantum photonic waveguides. 



PACS numbers: 42.50.-p, 42.70.Qs, 42.81.-i 



2 



1. Introduction 

The controlled generation of quantum entanglement is an essential process for 
performing a wide range of tasks in the field of quantum information |1]. 
Quantum communication protocols such as teleportation [2], key distribution [3] and 
dense coding |1] all require the generation of entangled states. Entanglement is also 
generated during quantum computation [5] and in the simulation of many-body quantum 
systems [6j. So far, small-sized entangled states have been generated in a wide range 
of physical setups, with photonic systems representing one of the most promising due 
to their speed and flexibility. Recently much attention has been focused on generating 
multiqubit entangled states such as cluster [TJ [8] , graph [9] and Dicke [TOj states, as well 
as those with more unusual structures and correlations [11] in a probabilistic fashion. 
Schemes that allow for the deterministic generation of smaller two-qubit entangled states 
have also been realized [121 ttSl HI] ? although low generation rates make these approaches 
challenging at present for building up to larger multiqubit entangled states. Photonic 
setups that generate states probabilistically, although inherently non-scalable, provide a 
readily available test-bed for probing the unique properties of quantum systems. Current 
experiments are however limited to entangled states of ten-qubits or less [15] . Improving 
photon generation rates using new types of sources and understanding better the 
practical requirements for generating high-quality entanglement between these sources 
may open up access to even larger entangled states with more complex structures. This 
would enable the testing of quantum protocols and probing physical phenomena that 
only more sizable quantum systems are able to support. 

In this work we report the first experiment to fuse photons from two independent 
photonic crystal fiber (PCF) sources into polarization entangled states. The fusion 
operation we demonstrate could be used to generate larger multiqubit entangled states, 
such as cluster and graph states for use in future quantum photonic technologies [T6] . 
including quantum communication and computation. We also introduce a novel method 
to characterize the fusion process that could be applied to a variety of different setups, 
including bulk, waveguide and on-chip photonic systems. While the fusion process is 
probabilistic in our experiment, we comment on how it may be made deterministic in 
the long term for the purposes of building large scalable entangled quantum systems. 
Our PCF sources have a very high photon generation rate compared to commonly 
used nonlinear crystal sources |l7j. Moreover, they make use of a zero-slope section of 
the phase-matching curve to produce photons in an intrinsically pure state, enabling 
good quality quantum interference without the need for spectral filtering (with inherent 
loss) resulting in high count rates [TTl HE]. We detect 90,000 coincidences per second 
from photon pairs produced with a pump laser power of 10.5mW. In this context, we 
investigate the fusion of photons from two independent sources in detail and provide a 
full quantum process reconstruction. We do this by developing a theoretical model 
to describe the entangling operation that includes the main causes of imperfection 
between the photons being fused. We find that our model fits well with the experimental 
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results. Our work shows how non-ideal settings limit the success of the fusion operation 
and provides detailed information about the practical requirements for using it to 
generate large multiqubit entangled states in bulk [171 HE] and on-chip quantum photonic 
settings [HI [201 ED [221 [22] . The fusion is performed with a fiber polarizing beamsplitter 
(FPBS), which is better suited than the bulk optics equivalent for integration into 
more complex schemes and scaling up to generate larger entangled states, where space 
constraints and coupling stability become important considerations. As it is waveguide- 
based, it is also interferometrically stable and well suited to the generation of resource 
states that are to be transmitted and received over fiber networks, for use in distributed 
quantum communication and networking protocols. 

2. Experimental setup 

Our setup for fusing photons from two independent PCF sources is depicted in Fig. [H 
A Ti:sapphire laser emits 8 nm pulses at a wavelength of 724 nm with a repetition rate 
of 80 MHz, which are then filtered to 1 nm. The pulses are split at a 50:50 beamsplitter 
(BS) and rotated to horizontal polarization by halfwave plates (HWPs) in both arms, 
then passed through polarizing beamsplitters (PBSs) and launched into the two PCFs. 
The fibers each produce a pair of photons (signal and idler) polarized orthogonally 
to the pump via four-wave mixing with nondegenerate wavelengths of 625 nm for the 
signals and 860 nm for the idlers. While these wavelengths can be tuned by changing the 
wavelength of the pump pulse, the photons are only emitted in an intrinsically pure state 
at this section of the phase-matching curve [TTl [18]. A 90° twist in the fiber ensures 
these photons exit the PCFs horizontally polarized, with an aspheric lens directing 
them into the polarizing beamplitters (PBSs), which transmit them. The PBS helps 
to filter out the pump pulses, which exit the fiber vertically polarized, and unpolarized 
background such as Raman scattering. An isolator (ISO) at the output of the laser, 
which allows only one-way transmission, blocks any reflected pulses from entering back 
into the laser. The photon pairs are then separated into different paths using dichroic 
mirrors (DM). The separated idler photons pass through a long- and short-pass filter 
tilted to give a tunable transmission window of 4 nm FWHM, which transmits the idler 
whilst removing as much Raman background as possible. They are then collected in 
single mode fibers followed by additional PBSs which remove background from reflected 
light in the PCFs. Multi-mode fibers are used to couple the idler to the detectors, where 
detections at D3 and D4 are used to herald the generation of the signal photons. The 
signal photons are rotated to diagonal polarization by HWPs and pass through a 40 nm 
FWHM bandpass filter. This bandwidth is large compared to the signal photon because 
the intrinsically pure state phase-matching makes narrow filtering unnecessary, and the 
filter only needs to remove any remaining light from the bright pump beam. They are 
then collected into single-mode fibers to guarantee optimal spatial overlap on the FPBS. 
Fiber polarization controllers (PCs) compensate for the effects on the polarization of 
strain-induced birefringence in the fibers. A quarter wave plate (QWP), HWP, QWP 
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Figure 1. (a): Experimental setup (b): Fusion operation used on two linear cluster 
states to make a new linear cluster state with an attached qubit (qubit 1) that can be 
used for subsequent fusions after a Hadamard operation is applied or for decoherence 
protection [23] . (c) : Fusion operation used on two linear cluster states to make a new 
two-dimensional cluster state. 

chain on mode 2' compensates any unwanted phase from the FPBS. Multi-mode fibers 
are used to couple the photons to the detectors. The detection of a photon at detector 
Dla (or Dlb) with detector D2a (or D2b) together with detections at detectors D3 and 
D4 are used in a fourfold coincidence circuit to register a successful fusion. Here, the 
use of two detectors on each signal mode enables active filtering of higher-order photon 
emissions from four-wave mixing in the PCFs. 

3. Fusion operation 

The fiber polarizing beamsplitter (FPBS) depicted in Fig. [1] transmits horizontal 
polarized photons (H) and reflects vertical polarized photons (V). Ideally it can be 
described by a unitary operation with respect to the reduced two-photon basis Tp : = 
{\HH) , \HV) , \VH) , \VV) , \{HV)Q) , |0(W))}, where \{HV)Q) {\Q{HV))) represents 
two photons of different polarization in the first (second) mode and the vacuum, 
|0) = |fac), in the second (first). In the case that we 'monitor' the number of photons 
present in the output ports, i.e. we gain information about the photon number but 
not the polarization, then the FPBS is described by the non-unitary trace-preserving 
channel £{p) = J2i=o -^iP-^l ^ with J2l=oM-^i — ^- Here the Kraus operator = 
\HH){HH\ + \ VV){VV\ + \HV){0{HV) \ + \ VH){{HV)0\ describes an operation into 
the 'coincidence basis', or subspace of the Hilbert space with states having a single 
photon in each output mode Cb := {\HH) , \HV) , \VH) , \VV)}, where Cb C Tp. The 
second Kraus operator Ei = \0{HV)){HV\ + \ {HV)0) {VH\ describes an operation out of 
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the subspace defined by the coincidence basis Cf, and into a subspace of the basis Tp with 
states having two photons in one of the output modes. As all operations are limited to 
the two-photon subspace defined by Tp, both Eq and Ei do not contain any basis states 
outside this subspace, such as \{HH)0). Furthermore, for input states restricted to the 
coincidence basis one can drop the last two terms from Eq and by monitoring the output 
of the FPBS in the coincidence basis only, allowing these states to be transmitted and 
rejecting all other cases, one can also drop the Ei operator to give the non-unitary non 
trace-preserving channel £f{p) = EopE^. This channel describes the combined action 
of the FPBS and monitoring in the coincidence basis with the Kraus operator 

Eo = \HH){HH\ + \VV){VV\. (1) 

This is a parity check operation [251 EE] in the sense that only states where both of the 
photons have the same polarization (even parity) are allowed to be transmitted. This 
operation can be used as a fundamental component in a variety of quantum protocols, 
such as quantum error correction [1], entanglement purification [26] and filtering [271128]. 
It also forms part of an efficient method to generate cluster state resources that can 
be used to carry out quantum computation [29j. Indeed, the parity check is used in 
Type-I fusion for efficiently building one-dimensional |30] and two-dimensional cluster 
states [30l [31] and can also be considered as a basic fusion operation; if successful, 
it matches Type-I fusion, but without the loss of any photons involved [H [32] . In our 
experiment, however, in order to 'monitor' the output photon number in the coincidence 
basis we detect the photons by destroying them. Thus our parity check, or fusion 
operation, is a postselected version and not scalable in its current form, unlike Type- 
I fusion |30j. This is due to an exponential decrease in the total success probability 
when using many of these operations. It should be noted that active monitoring in 
the coincidence basis without destroying the photons can be achieved in principle using 
a non-demolition type measurement of the photon number [33], making this approach 
scalable using the techniques of Ref. pO]. However, such measurements require the 
photons to pass through media with a large nonlinear response at the single-photon 
level, which is technically challenging at present and beyond the scope of this work |34j . 

Ideally, the FPBS channel Sp from Eq. ([T]) takes two photons in the state 
Pin = \in){in\, with \in) = |+) |+) and |+) = + |V^)), and fuses them into 

the maximally entangled state = {l/\/2)[\HH) + \VV)], equivalent to a two- 

qubit cluster state under a local Hadamard rotation on the first qubit. The above 
fusion operation occurs with success probability po = TT^EopinE^) = 1/2 due to the non 
trace-preserving nature of the channel. Note that regardless of the success probability 
of a given input state pm being transmitted into the coincidence basis, the channel Sp 
always acts with unit probability. Either a coincidence is measured and the input state 
is transmitted into the coincidence basis, or no coincidence is measured and the input 
state is not transmitted into the coincidence basis. Moreover, if both photons represent 
the end qubits of two seperate linear cluster states, as shown in Fig. [1] (b) , then after 
passing through the channel Sp the cluster states are fused together; here the two linear 
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cluster states involved are joined with one of the photons representing the middle qubit 
between the edges bonding to the two clusters and the other photon representing a qubit 
directly attached to it |9]. If failure occurs one tries again with the next qubits in both 
clusters [30]. Two-dimensional cluster states can also be formed using the channel £p as 
shown in Fig. [1] (c) , or in a scalable way based on the techniques of Refs. |30] and |31] . 

4. Fusion interference 

In order to check that our FPBS is implementing the correct fusion operation Ep and 
coherently interfering two input photons in the required spatio-temporal manner we 
set the input state of the two signal photons to \in) = |+) |+) and modify a time 
delay 6t between the input photon pulses using a translation stage, as shown in Fig. [H 
This allows us to maximize the coincidence probability temporally, with the FPBS 
single-mode input fibers ensuring good spatial overlap. At the detectors we measure 
the coincidence probability for the |+) |+) state population, which for the expected 
entangled state would give 1/2. However, the state is produced from the 
fusion process acting on the input state |m) with a probability 1/2, so we expect the 
total coincidence probability to be 1/4. On the other hand, for an incoherent interference 
of the two photons, i. e. when St is larger than the mutual coherence time of the incoming 
photons Tc, but less than our detector coincidence time- window Tcoinc, the probability 
for a coincidence is 1/8 upon considering all possible outcomes of the non-interfering 
photons. Thus, as we modify the time delay 6t in our setup, we look for an 'antidip' 
or peak in the coincidence probability (number of coincidences). This is a Hong-Ou- 
Mandel type interference effect [35] and a rigorous time-dependent derivation of the 
quantum interference phenomenon leading to this antidip in the coincidence probability 
is outlined in Appendix A. Here one finds that the probability for a coincidence at a set 
time delay of 6t is given by 

Pcoinc(5r) = i(e-(^-/-)' + l), (2) 
o 

where at is the pulse duration of the signal photon, defined as 2v21ii2/Aw, and Acu 
is the full width at half maximum of the signal's spectral intensity. In our experiment 
at — Ips and Au corresponds to AA ~ 0.5nm. For complete temporal overlap 5t = 
and for coherent interference we expect a coincidence probability of 1/4. On the other 
hand for large 6t, \6t\ ^ Tc, where ^ at, and with the coincidence time- window larger 
Tcoinc ^ |5t| we have incoherent interference and the coincidence probabihty drops to 
1/8 as mentioned above. 

In Fig. [2] (a) we show the expected number of coincidences A^'coinc('5T) obtained from 
Eq. (I2l). Here, NcomciSr) = NavPcoinci^T), where we have multiphed Pcoinc(5T) by the 
average number of counts in our experiment when \6t\ ^ in order to scale it correctly. 
We have A^'av = 401. The function A^coinc(^T) is plotted as the solid upper curve. In 
the lower dotted curve we include a possible frequency mismatch between the central 
frequencies of the two signal photons, and cu^, corresponding to AAm = 0.06 nm 
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Figure 2. Coincidences and tomographic reconstruction of tlic density matrix for the 
fused state resulting from input state |+) |+). (a): Fourfold coincidences measured 
as a function of the delay time 5t (divided by the signal pulse duration at) between 
the heralded signal photons clearly showing the expected antidip around zero delay. 
Using a relatively low pump power of 5.3 mW per fiber in this case, we obtain 23,000 
pair coincidences/s per fiber and a maximum of 1.25 fourfold coincidences/s. The 
top curve is the ideal theoretical dependence and the dashed line includes a possible 
frequency mismatch corresponding to AAm = 0.06 nm for the signal photons. The 
lower solid line is a best fit of the data to the function A'av(poe~'''^^^'^'^ + 1)/8, finding 
Po = 0.61 and A'av is the average number of counts in our experiment when St ^ at- 
(b): Experimental real part of density matrix, (c): Experimental imaginary part, 
(d): Ideal real part for (e): Ideal imaginary part for 

coming from the spectral resolution of our spectrometer. This is calculated by carrying 
out the necessary integration over the coincidence time-window for a spectral-temporal 
dependent version of Eq. ([2]) derived in Appendix A as Eq. ( lA.SP . 

The data points from our experiment are also shown in Fig. [2] (a), along with 
a best fit to the curve defined by N^oind^'^) = ^a-vf'IoLci^''')^ where PcoLd^'^) = 
{poe~^^'^^'^*^ +l)/8, withpo = 0.61, giving an antidip visibility of 61%. Here, the visibility 
is defined as the percentage ratio of (A^max — A'min)/A^min- In our experiment we have 
Tcoinc = 3 ns, Trep = 12.5 US and 6t is varied using a translation stage on the pump beam 
before entering one of the PCFs, as shown in Fig. [1] (a). The errors are calculated using 
a Monte Carlo procedure with Poissonian fluctations in the count statistics [37]. The 
main causes of deviation from the ideal case are the spectral mixedness of the signal 
photons from our sources [18] , leading to broadening of the coincidence profile [36] and 
spatial mismatch between the interfering signal photons at the FPBS. We discuss later 
how these effects can be incorporated into our theoretical model. Note that photon loss 
in general affects the overall number of coincidences regardless of the time delay St, 
leaving the visibility of the antidip unaffected. 

In Fig. [2] (b) and (c) we show a tomographic reconstruction of the experimentally 
fused entangled state of the two signal photons at 5r = 0. This is obtained by measuring 
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Table 1. Properties of the fused state as the pump power is increased. 



power (mW) 


4-fold rate (s^i) 


pairs/pulse 




C 


P 


5.3 


3.2 


0.037 


0.740 ±0.007 


0.550 ±0.014 


0.63 ±0.01 


7.9 


9.8 


0.064 


0.677 ±0.006 


0.392 ±0.012 


0.520 ±0.007 


10.5 


36.4 


0.103 


0.606 ±0.007 


0.265 ±0.015 


0.448 ±0.008 


13.2 


77.8 


0.160 


0.554 ±0.006 


0.15 ±0.01 


0.392 ±0.005 


14.8 


111.6 


0.193 


0.520 ±0.004 


0.07 ±0.01 


0.359 ±0.003 



the qubits in all combinations of the bases {\H/V) , |+/— ) , \R/L)} [37]. Ideally this 
state is the maximally entangled state l^"*") shown in Fig. [2] (d) and (e). At 3.2 four- 
fold coincidences per second our experimental state has a fidelity of = 0.74 ± 0.01 
with respect to l^"*"), entanglement quantified by the concurrence of C = 0.55 ± 0.01 
and a purity of P = 0.63 ± 0.01. Some background can be seen to be present in the 
state from the non-zero \HV) and \VH) terms in Fig. [2] (b), which we attribute to 
higher-order photon emission, other background processes such as Raman scattering, 
and imperfect polarization operations in our measurements. The background also gives 
rise to several non-zero terms in Fig. [2] (c). In Appendix B we analyse the effects of 
higher-order photon emission from the four-wave mixing process and find that at this 
count rate we expect the fidelity to be limited to F^+ = 0.80. By increasing the pump 
power to 14.8 mW we achieve a four-fold coincidence rate of 111.6 per second with 
a fidelity of = 0.520 ± 0.004, a non-zero concurrence of C = 0.07 ± 0.01 and a 
purity of P = 0.359 ± 0.003. In Table [1] we show how various parameters change as 
the pump power is increased and higher-order emission becomes increasingly significant. 
These generation rates and corresponding state properties show the potential of using 
PCFs for the controlled generation of larger multi-qubit entangled states for carrying out 
quantum protocols and investigating quantum phenomena with desired state qualities 
and success rates. 

Spatial-temporal mode mismatch, loss and higher-order photon emission are sources 
of error leading to deviation of our experimental data from the ideal case in both the 
coincidence probability and state tomography. While an explicit analysis of all factors 
leading to these effects is beyond the scope of the current work, we will show how it 
is possible to develop a theoretical model to describe the non-ideal fusion operation 
generating entangled states and connect it to the experimental antidip coincidence 
probability curve shown in Fig. [2] (a). 

5. Fusion process tomography 

Ideally our fusion operation should act as the non trace-preserving channel £f{p) = 
EopE^, with Eq defined in Eq. (Q. However, as the experimental fusion is not ideal, 
we must find a more approriate channel description. Assuming polarization changes 
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Outpii) States Output States Output States 



Figure 3. Input and output states in tlie probing of the fusion operation, (a): 
Experimental coincidences for the Z ^ Z basis, (b): Ideal probability for the Z ^ Z 
basis. The fidelity Fz^z is given by the ratio of correctly transmitted photon pairs 
to total number of transmitted pairs by the fusion operation for input basis Z and 
output basis Z. Similarly for Fx~^x and Fx^y- (c): Experimental coincidences for 
the X —i' X basis, (d): Ideal probability for the X —i' X basis, (e): Experimental 
coincidences for the X Y basis, (f): Ideal probability for the X — F basis. In all 
plots, each of the four-fold coincidence points was taken in 621s. 



from input to output states are negligible (see Fig. [3] (a)) we can write the experimental 
fusion operation as [28] 

ST'iP) = Y.<mEnpEl (3) 

n,m, 

where the operators En are given as 

Eo =\HH){HH\ + \VV){VV\ 
E,, = \HH){HH\-\VV){VV\ 
E,y = \HV){HV\ + \VH){VH\ 

E,, = \HV){HV\-\VH){VH\, (4) 

which form an orthogonal operator set with respect to the Hilbert-Schmidt inner 
product. £"0 is the ideal operation which occurs with probability Xo,oy the operator E^z 
describes a phase flip error with probability Xzz zz^ the operator E^y describes leakage 
of odd-parity states {\HV) and \VH)) through the fusion device with probability x. 



and finally the operator E^^ describes leakage and a phase flip with probability xf 



^xy,xy 
F 

^xx,xx' 

Thus we have n,m E {0, zz,xy,xx} and J2nXn,n = 1- We can also define the basis 
fidelity [38] for input basis states i and output basis states j as 
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where for Fz-,z we have {\HH) , \HV) , \VH) , \VV)} {\HH) , \VV)}, for Fx^x we 

have {|++) , I )} ^ {|++) , I )} and {|+-) , |-+)} ^ {|+-) , |-+)}, for Fx^y we 

have {|++) , I )} ^ {\LR) , \RL)} and {|+-) , |-+)} ^ {\LL) , \RR)}. Substituting 

the above i and j basis states into Eq. ([5]) and using the orthogonahty of the En operators 
appearing in S"^^ of Eq. ([3]), the basis fidehties are found to be equivalent to the elements 
of the X matrix as follows [28] 

Fz^Z = Xo,0 + Xzz,zz 
Fx^X = Xo,0 + Xxx,xx 

Fx^Y = Xo,0 + Xxy,xy (6) 

From these we then have the relation for the process fidelity of the fusion operation, 
F^^ = Xoo describing the probability for successful fusion as 

Fp' = \{Fz^z + Fx^x + Fx^Y - 1), (7) 

where we have used J2n Xn,n — 1- The basis fidelities can be found from our experiment 
using the relation 

Tout 

1& LliC liUlllUBl Ui tllllCO LiiC uuipuL sua 

rout 



as described in Appendix C. Here, N°^^j^ is the number of times the output state is 



measured to be \je) when input state \ik) is sent through device and Tj^. = Xljp^ip.ifc- 
Eq. (IH]) essentially means we can define the basis fidelities Fi^j to be the ratio of the 
total number of transmitted output states in the correct basis to the total number 
of transmitted states. In Fig. [3] we show the experimentally measured number of 
transmitted states (coincidences) for the input and output basis defined by the fidelities 
Fz^z, Fx~^x and Fx^y- We find Fz^z = 0.958 ± 0.001 from the data shown in 
Fig. |3] (a), with the ideal transmission coincidence probabilities for Z Z shown in 
Fig. [3] (b). From the data corresponding to Fig. |3] (c) we find Fx^x = 0.768 ± 0.001. 
The ideal transmission coincidence probabilities for X — )■ X are shown in Fig. |3] (d). 
Finally, from the data corresponding to Fig. |3] (e) we find Fx^y = 0.759 ± 0.001, with 
the ideal transmission coincidence probabilities for X Y shown in Fig. |3] (f). From 
these three fidelities we then find, using Eq. ([7]), a process fidelity of our fusion operation 
of Fp^ = 0.743 ± 0.001. 

The entanglement capability of our fusion operation can be defined as the maximum 
entanglement that can be generated from it using input product states. We can quantify 
this in terms of the concurrence C and we have the lower bound for the entanglement 
capability of the fusion operation, quantified by the concurrence as C^; > 2Fp^^ — 1, as 
shown in Appendix D. Based on the above process fidelity of our experiment we have 
Ce > 0.485 ± 0.002. Note that both Fp^^ and Ce are consistent with the fidelity and 
concurrence calculated directly from the reconstructed density matrix of the last section 
for the rate 6.62 events per second at 5r = 0. A small deviation is due to the assumption 
of no polarization change from input to output states in the model used to describe the 
experimental fusion process S*^^. 



11 



6. Non-ideal temporal setting 

So far we have considered our fusion operation at the maximum visibihty 6t = 0. In the 
ideal case oi 5t — we have that the state |+)-^ |+)2 is transformed into the entangled 
state \4>'^)ii2i with probability 1/2. However, for \6t\ 3> Tc we obtain the incoherent 
mixture p = \{\HH){HH\ + \VV){VV\)v2' with probability 1/2. Here, the photon in 
mode V is detected before the photon in mode 2' (or vice-versa) and even though the 
time period \5t\ is outside the coherence time of the individual photons Tc, the larger 
coincidence window Tcoinc means we still register the state as being a valid output state 
from the fusion process. This loss of coherence in the output state suggests the process 
might be described by a phase damping channel acting on the output photons. This 
can be incorporated into the fusion channel by writing the total action of the fusion as 
^riPin) = £pD{£F{Pin))) where £pd is a phase damping channel, the form of which we 
derive next. The same channel might also describe other forms of decoherence, such as 
that resulting from imperfect spatial-mode overlap or spectral distinguishability between 
the interfering photons. These parameters are challenging to set in an experiment so it 
is useful to understand how the fusion operation performs under such conditions. By 
developing a model that describes the effect of non-ideal control we may then quantify 
the requirements for a given level of performance of the fusion. 

For a single qubit, a phase damping channel with an arbitrary time dependent 
dephasing function /(5r) G [0, 1] can be described by the channel £{p) = KipKj , 
with the Kraus operators 

which satisfy K^Kq + kIKi — 1. Usually the function /(5t) = e"^^'^ is considered, 
where F represents some decay rate and St is a positive exposure time. Here we wish to 
determine the form of /(5r) for our fusion operation initially using St as the temporal 
parameter. Considering two qubits subject to such a type of arbitrary phase damping, 
we have the overall channel £pd{p) = Y^i=o^iP^li where 

ki = apt <S> 
k2 = apa^ (8) 1 

ks = (8) a„ (9) 
with a = [(1 + /(5t))/2]V2 ^ ^ [(i _ j(5^))/2]V2. ^^^^^ ^^^^^ 

channel St{p) — £pd{£f{p)) describing the non- ideal fusion and phase damping due to 
(5t 7^ in the {E^} basis as 

^t(p) = xi^k.KpElkl = J2 xlrrApEl (10) 

i,n,m n,m 
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The second equality is possible because H and do not change the polarization, so 
that we can rewrite the combined operators more conveniently in terms of only £"0, E^z, 
Exx-, Exy, and the elements of a new matrix using the identities H ® H = £'0 + Exy, 
H ® = Ezz — Exx, cr^ ® 11 = Ezz + Exx and az ^ (Tz = Eq — Exy. In particular, we 
have that the total success probability for the fusion operation is given by Fp = Xo — 
(a^ + (3'^)Xoo + (2tt^/3^)xfz zz- Substituting in for a and (3, with /(5r) arbitrary, we have 

= lixlo + XL,J + \f{^rf{xlo - Xlzj- (11) 
Similarly, we have 

T ^ ( F ^ F \ 1 £( ^ \2 { F F \ 

Xzz,zz — 2^^0,0 + Xzz,zz) ~ 2''^ ' 1X0,0 ~ Xzz,zz) 

\^ =-ix^ ) + -f(5Tf(\^ -Y^ ) 

A.xx,xx 2\A.xx,xx ' A.xy,xyJ ' \" J \A,xx,xx A.xy,xyJ 

X^ =-(x^ +X^ )--fi5T)HY^ -Y^ ). 

'^xy,xy C2\^xx;xx ' ^xy^xyl c^'^ \ I \^xx,xx ^xy,xyJ 

(12) 

We can find the expected form of the phase damping function fiSr) by considering the 
coincidence probability Pcoinc given in Eq. ([2]) for the case of ideal fusion. Here, the 
fusion operation acts on the state |+) |+) and produces the state 

p = ^[\HH){HH\ + \HH){VV\ + \VV){HH\ + \VV){VV\l (13) 

with probability 1/2. After phase damping we then end up with the state p' = £pd{p) 
given by 

p' = ^[\HH){HH\ + f{5Tf\HH)iyV\ + f {5Tf\VV) {H H\ + \VV){VV\i (14) 

Measuring the population of p' in the state |+) |+) gives (1 + /(5r)^)/4. By including 
the factor of 1/2 due to the probabilistic nature of the fusion producing the state in 
Eq. <^ we have that Pcoinc = (1 + /((5r)2)/8. Using Eq. <^ we are able to then 
make the correspondence /(5t) = e"^'^'^/'^*^^/^. In Fig. H] (a) we show the expected 
dependence of the total process fidelity Fj, and entanglement capability with the time 
delay 5t in the case of our experimental fusion at 5t = 0. The data give a close match 
to the expected dependence. For definiteness we choose positive St with the temporal 
dependence being symmetric. In both of these we have used the phase damping function 

In Fig. m (b) we show the matrix elements for the total process of fusion 
and phase damping, as defined in Eqs. (fTTj) and ( IT2|) . The solid lines correspond to 
the expected dependence of the ideal fusion with phase damping and the dotted lines 
correspond to the expected dependence of our experimental fusion at 5r = with 
phase damping applied. In both, the phase damping function is /(5r) = e"^^'^''"'*-'^/^. 
Again, the data closely match the /(5t) model. Note that for ideal fusion, the elements 
X^x,xx X^y^xy ^^^o always, regardless of the amount of phase damping. However, 
Fig. m (b) shows that this is not the case for nonideal fusion, as in our experiment. Note 
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Sr/<ri 



Figure 4. The total process fidelity, entanglement capability and elements of the 
process matrix, (a): The total process fidelity Fp and entanglement capability Ce 
dependence on the time delay 5t. The solid (dotted) line shows the decay expected 
from the theoretical phase damping model with the experimental fusion values taken 
at zero time delay for Fp {Ce)- The experimental values are plotted as triangles 
(circles), with errors negligible on the scale of the graph, (b): The elements of the 
total process matrix as they depend on the time delay St. The upper (lower) solid 
line corresponds to the element Xo.o ix^z zz) from the theoretical model with perfect 
fusion. The upper (lower) dashed line corresponds to the element Xo ^ ix^z zz) from 
the experimental data dX 5t — Q extrapolated using the function /{St). The lower 
lines show the elements Xxx.xx and Xxy.xy Fo'" ideal fusion these are zero regardless of 
the time delay. The experimental values are plotted as circles, triangles, squares and 
diamonds for 00, zz, xx and yy respectively, with errors negligible on the scale of the 
graph. Note that all diagonal elements sum to one regardless of the time delay 5t. 

also from Fig. H] (b) that in the case of our experimental fusion at 5t = 0, the effects of 
imperfections such as spectral mixedness and spatial mismatch can be almost entirely 
described in terms of an ideal fusion operation (solid lines) that has been phase damped. 
Indeed, at Sr/at ~ 1 in Fig. H] (b), one can see that the phase-damped ideal fusion 
matches the non-damped experimental fusion at 5r = 0. Thus to a good approximation 
iXxx,xx Xxy,xy negligible), we can model our experimental fusion operation at 5r = 
as being a phase damped ideal fusion operation. For (5r > an additional phase damped 
channel with the function /{St) is then applied to each qubit. 

7. Summary 

In this work we reported an experimental demonstration of the fusion of photons from 
two independent photonic crystal fiber sources into polarization entangled states. We 
introduced and carried out a novel method to characterize the fusion operation via 
quantum process reconstruction. To do this we developed a theoretical model to describe 
the entangling process which included the main imperfections. We then showed that 
our model fitted well with the experimental results. This work highlights the need for 
accurate control of spatial and temporal properties for the success of photonic fusion 
and provides detailed information about the practical requirements and limitations of an 
operation that may be used in emerging quantum photonic technologies. Future work 
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will be to optimize our fusion operation and use it to generate larger high-quality multi- 
qubit entangled states, including the use of additional degrees of freedom such as path 
and frequency encoding, for carrying out computations and communication protocols. 
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Appendix A: Derivation of fusion interference antidip 

Starting from the initial state 1-^^)2 ^^e signal photons in modes 1 and 2, 
produced from the PCF sources as shown in Fig. [T|, and finishing with the measurement 
of the coincidence probability for the state \H)^^ \H)2a detectors Dla and D2a, 

we have the probability density for a photon to be detected in mode la at time to and 
another in mode 2a at time 1^ + t given by 



Pcoinc(to,r) = ,2(ra|^Hi„(^o)^H,„(io + r)i+,Jto + r)i+^^(to)|m)i2 (A.l) 



where the scalar photon-unit time-dependent electric field operators are defined as 



with a\{uj) and dj(u;) the photon creation and annihilation operators respectively, which 



where (piloJi) is the spectral amplitude of a single photon pulse in mode i, normalized 
so that / dijj(f)*{ijJi)(j)i{ui) = 1. The initial state \H)^ \H)2 is rotated into the required 
input state \in) = \+)2 of the FPBS by HWPs that produce the transformations 



a^H^itui) {l/V2)[a^H^iuji) + a\,^iui)] and a^H^iu2) ^ {1/ V2)[a^jj^{cu2) + al^^M]- The 



action of the FPBS, combined with the HWP-QWP-HWP chain on mode 2 (for phase 
correction), produces the transformations 




(A.2) 




al^{u;2) -m{/^,(u;2) 



(A.4) 



Next, the HWP and QWP on each mode produce the transformations 
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The PBSs that follow produce the transformations 

Thus, only horizontally polarized photons will be detected at detectors Dla and D2a. 
These last operations (the HWP, QWP and PBS) together with the detections in modes 
la and 2a are equivalent to measuring the coincidence probability for the state 1+)-,^, |+)2/ 
in the output of the fusion. Therefore Pcoinc(^05 t) in Eq. (lA.ll) represents the coincidence 
probability density for the state |+)-^, |+)2/- 

Substituting all the above transforms into Eq. ( lA.Sp . then substituting this into 



Eq. ( ]A.1|) with the definitions given in Eq. (1A.2P and carrying out the necessary 



integrations [36j one finds 

Pcoinc(to,r) = ^\Ci{to + r)C2(to) + Ci(to)C2(to + (A.5) 
lo 

where Q{t) = J dui(f)i{cOi)e~^'^'^^ is the spatio-temporal modefunction for mode i. 
Choosing Gaussian single-photon pulses for each of the input modes, with a time delay 
of St between the peaks of the pulses, we have 

Cl(t) = (2/7r)l/^e-(*-^"/2)2-ia;?t ^j^Q^ 

(^(t) = (2/vr)i/4e-(*+^"/2)2-i-2°t. (A.7) 

Here, is the central carrier frequency for mode i expressed in units of l/crt and the 
times t and 6t are expressed in units of at {at is the pulse duration of the signal photon, 
defined as 2 a/2 In2/Aa;, and Auj is the full width at half maximum of the signal's spectral 
intensity). Substituting these expressions into Eq. f lA.5p and integrating over all possible 
detection time to, i.e. dto, one finds the probability density for a coincidence at a 
time duration r with a set delay of 5t given by 

Pcoinc(r, 5t) = _ (cos[r(u;° - a;°)] + cosh[25r r]). (A.8) 
V647r 

Setting the central carrier frequencies the same for the moment, uj^ = U2, we then 
integrate Eq. (lA.SP with respect to r over the coincidence window. Here, the integration 
is r""'""^^ allowing for either signal photon to be detected first. For Trep > Tcoinc ^ 

''"coinc/ ^ " 

|5r|, where r^ep is the time between the pump pulses (and therefore possible signal 
photons produced from four- wave mixing in the PCFs), one finds the coincidence 
probability 

Pcoinc(5r) = l(e-(^-/-)' + l). (A.9) 
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Where we have now included at exphcitly. Although we have neglected a time-dependent 
evaluation of the idler photons in the above analysis, they can also be included. Indeed, 
it is straightforward to check that such a calculation does not change the result of 
Eq. (]A.9|) . given that the idler photons are assumed to be in a product state with the 
signal photons. 



Appendix B: Higher-order emission analysis 

The state generated by four-wave mixing in one source can be written as [39] 

= ATdO, 0),^, + a \H, H)^^^ + a' \2H, 2H)^^^ + ©(a^)), (B.l) 

where A/" is a normalisation constant and |ap = n is the mean number of signal-idler 
pairs generated in a pulse. When the idler detector registers a click, the heralded density 
matrix of the signal mode becomes 

where the action of the detector is described by the positive operator-valued measure 
{Ilciick, Ilno-ciick = H - Ilciick}, with Ilciick = Y.n=o(^ " (1 - vT) l^H) ^ {nH\ and 77 is the 
lumped detector efficiency of registering a click given a single photon is input into the 
mode that it monitors. From Eq. ( 1B.2|) we have 

Ps = \H) {H\ + n7]2 \2H) {2H\ + 0{n^)), (B.3) 

where = 1 — (1 — 77)". We now write the heralded state of the two sources to first-order 
in n as 

pu = Af{\H, if) ,2 {H, H\ + 2nj \H, 2H),^ {H, 2H\ + 2nj \2H, H),^ {2H, H\), (B.4) 

where 7 = ri2/2rii. The signal photons from pi2 are first rotated from horizontal to 
diagonal polarization then input to the fiber polarizing beamsplitter (FPBS). Using the 
transformations given in Eq. flA.4p of Appendix A we can write the output from the 



FPBS by making the following substitution in Eq. flB.4p 

Where we have now switched from the Schrodinger to the Heisenberg picture. By 
removing the terms which do not lead to a possible coincidence of detection clicks 
between the two signal modes, then performing the polarization rotations used in 
analysing each measurement basis {\H/V) , |+/— ) , \R/L)} and switching back into the 
Schrodinger picture to determine the contribution to the coincidence counts of the 0{n) 
terms in Eq. (lB.4p . we can evaluate the effects of higher-order emissions in our data. 
For instance, we find the visibility of the anti-dip in Fig. O up to 0(n), is limited to 

^ 1 - 4717 + 6n7^ 
1 + 6rl7 + 37172 



17 



In addition, for the lowest power in Table 1 where n = 0.037, taking 77 = 0.1 and 
considering up to 0(ji) terms, the expected fidelity of the output state with respect to 
10^) is F^+ = 0.80. This assumes the fusion process is ideal, and so represents an upper 
bound. 



Appendix C: Derivation of basis fidelities 

The basis fidelities can be found from our experiment by noting the relation 

{MSr'{\^k){^k\)m = ^, (C.l) 

where N-^ is the number of times input state \ik) is sent through the device and N°^j^ 
is the number of times the output state is measured to be when input state \ik) is 
sent through device. We then have that 

1 A^°"* 



m 



Keeping the value of N-^ constant over the experiment for all possible input states 
\ik), i.e. Nl^ = A^*", and noting that Tr(£^^^^(]l/4)) = 1/2 for any experimental process 
matrix x appearing in Eq. ([3]), we have that A^*" = where Tj^. = ^iyl^- 

In other words, for a given input basis the total number of transmitted states over 
that basis is equal to half the total number of input states (in this case 2iV*"). Therefore 
we have 

^i-^j = T — ■ ^ ' 

Appendix D: Relation of concurrence to fusion operation fidelity 

The concurrence, C, for a two qubit state pyi can be written as [ID] 

C(pi2) = max{0,- min Tr((|A)(ArVi2)}, (D.l) 

AeSl(2,C) 

where 1^4) denotes the unnormalized state {A® 1) |/), with |/) = and A is 

any matrix with det(A) = 1. In Eq. (ID.ip any choice for the state \Al) provides 
a lower bound on the concurrence, thus choosing \A) = y/2\ijj^), we have that 
{\A) {A\)^^ = 1 - 2|0+)(0+|, and we can rewrite Eq. flD?!]) as 

C(pi2) >max{0,2F<^+-l}, (D.2) 

where = Tr(|0+) (0+|pi2) is the fidelity of the state pu with respect to If 
we consider the state pi2 to be an output state from a channel S acting on two input 
product states, then C(pi2) provides a lower bound on the entanglement capacity of 
that channel. Choosing the input product state to be pm = |+)(+| ® |+)(+| and the 
channel to be our experimental fusion operation we have that pi2 = £l^^{pin) = 
Fp^^'l^"'") (0+1 + (1 — Fp^^)p'. Thus, = Fp^^ and we have the lower bound for 
the entanglement capability of the fusion operation, quantified by the concurrence as 
Ce > 2F;"^ - 1. 



18 



References 

[1] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge 

University Press, Cambridge (2000). 
[2] C. H. Bennett et ai, Phys. Rev. Lett. 70, 1895 (1993). 
[3] A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

[4] C. H. Bennett and S. J. Weisner, Phys. Rev. Lett. 69, 28812884 (1992); D. Brufi et at, Phys. Rev. 
Lett. 93, 210501 (2004). 

[5] D. Deutsch and R. Jozsa, Proc. R. Soc. Lond. A 439, 553 (1992); P. Shor, SIAM J. Comput. 26, 

1484 (1997); L. K. Grover, Phys. Rev. Lett. 79, 325 (1997); D. Colhns, K. W. Kim and W. C. 

Holton, Phys. Rev. A 58, R1633 (1998). 
[6] S. Lloyd, Science 273 1073 (1996); L Buluta and F. Nori, Science 326 5949 (2009). 
[7] P. Walther et at, Nature (London) 434, 169 (2005); N. Kiesel et ai, Phys. Rev. Lett. 95, 210502 

(2005); K. Chen et ai, ibid. 99, 120503 (2007); G. Vallone et ai, ibid. 98, 180502 (2007); Y. 

Tokunaga et ai, ibid. 100, 210501 (2008); R. Ceccarelh et ai, ibid. 103, 160401 (2009); D. N. 

Biggerstaff et ai, ibid. 103, 240504 (2009). 
[8] A.-N. Zhang et ai, Phys Rev A 73, 022330 (2006). 

[9] C.-Y. Lu et ai, Nature Physics 3, 91 (2007); S. M. Lee et ai, arXiv: 1105.5211 (2011). 
[10] N. Kiesel et ai, Phys. Rev. Lett. 98, 063604 (2007); R. Prevedel et at, ibid. 103, 020503 (2009); 

W. Wieczorek et ai, ibid. 103, 020504 (2009). 
[11] M. Radmark, M. Zukowski and M. Bourennane, Phys. Rev. Lett. 103, 150501 (2009); C.-Y. Lu et 

ai, ibid. 102, 030502 (2009); J. K. Pachos et at, New. J. Phys. 11, 083010 (2009); J. Lavoie et 

ai, Phys. Rev. Lett. 105, 130501 (2010); A. Chiuri et ai, arXiv: 1105.2109 (2011); L. Aolita et 

ai, arXiv: 1105.3598 (2011). 
[12] C. Wagenknecht et ai, Nature Photonics 4, 549 (2010). 
[13] S. Barz et at, Nature Photonics 4, 553 (2010). 
[14] H. Hiibel et ai, Nature 466, 601 (2010). 

[15] W.-B. Gao et at, Nature Physics 6, 331 (2010); X.-C Yao et ai, arXiv: 1105.6318 (2011). 

[16] J. L. O'Brien, A. Furusawa and J. Vuckovic, Nature Photonics 3, 687 (2009). 

[17] M. Haider et ai, Optics Express 17, 4670 (2009). 

[18] A. Clark et al, arXiv;1102.4415 (2011). 

[19] A. Politi et ai, Science 320, 646 (2008). 

[20] L. Sansoni et ai, Phys. Rev. Lett. 105, 200503 (2010). 

[21] M. S. Tame and M. S. Kim, Phys. Rev. A 82, 030305 (2010). 

[22] Q. Wang, G. Farrell and Y. Semenova, IEEE Sel. Top. Quant. Elec. 12, 1349 (2006); Q. Wang 

and G. Farrell, J. Opt. Soc. Am. A 24, 3303 (2007). 
[23] A. Crespi et al, arXiv:1105.1454 (2011). 

[24] R. Prevedel et ai, Phys. Rev. Lett. 99, 250503 (2007); M. S. Tame, M. Paternostro and M. S. 

Kim, New. J. Phys. 9, 201 (2007). 
[25] T. B. Pittman, B. C. Jacobs and J. D. Franson, Phys. Rev. A 64, 062311 (2001). 
[26] J.-W. Pan, C. Simon, C. Brukner and A. Zeilinger, Nature 410, 1067 (2001); J.-W. Pan, S. 

Gasparoni, R. Ursin, G. Weihs and A. Zeilinger, Nature 423, 417 (2003). 
[27] H. F. Hofmann and S. Takeuchi, Phys. Rev. Lett. 88, 147901 (2002). 
[28] R. Okamoto et al., Science 323, 483 (2009). 

[29] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 5188 (2001); R. Raussendorf, D. E. Browne, 

and H. J. Briegel, Phys. Rev. A 68, 022312 (2003). 
[30] D.E. Browne and T. Rudolph, Phys. Rev. Lett. 95, 010501 (2005). 
[31] Y. S. Weinstein, arXiv: 1 102:1916 (2011). 

[32] T. P. Bodiya and L.-M. Duan, Phys. Rev. Lett. 97, 143601 (2006). 

[33] N. Imoto, H. A. Haus and Y. Yamamoto, Phys. Rev A 32, 2287 (1985); P. Yeh, J. Opt. Soc. Am. 
B 3, 747 (1986); H-A. Bachor et ai, Phys. Rev. A 38, 180 (1988). 



19 



[34] C. Guerlin et ai, Nature 448, 889 (2007). 

[35] C. K. Hong, Z. Y. Ou, and L. Mandel, Phys. Rev. Lett. 59, 2044 (1987); R. Ghosh, C. K. Hong, 
Z. Y. Ou and L. Mandel, Phys. Rev. A 34, 3962 (1986); R. Ghosh and L. Mandel, Phys. Rev. 
Lett. 59, 1903 (1987). 

[36] T. Legero et at, App. Phys. B 77, 797 (2003). 

[37] D. F. V. James et al., Phys. Rev. A 64, 052312 (2001). 

[38] H. F. Hofmann, Phys. Rev. Lett. 94, 160504 (2005). 

[39] J. Fulconis. O. Alibart, W. J. Wadsworth and J. G. Rarity, New J. Phys. 9, 276 (2007). 
[40] F. Verstraete, Ph.D. thesis Katholieke Universiteit Leuven, 2002 (unpublished). 



